Universality of Uhrig dynamical decoupling for suppressing qubit pure dephasing and relaxation 
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The optimal A'-pulse dynamical decoupling discovered by Uhrig for a spin-boson model [Phys. Rev. Lett, 
98, 100504 (2007)] is proved to be universal in suppressing to 0(t'^*'\ the pure dephasing or the longitudinal 
relaxation of a qubit (or spin- 1/2) coupled to a generic bath in a short-time evolution of duration T. It is also 
found that for the purpose of suppressing the longitudinal relaxation, an ideal Uhrig ;r-pulse sequence can be 
generalized to a sequence consisting of the ideal one superimposed with finite-duration pulses satisfying certain 
symmetry requirements. 
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Introduction - A central topic in spin resonance spec- 
troscopy yj] is the decoherence of spins due to coupling to en- 
vironments, including the longitudinal relaxation of the pop- 
ulation and the transverse relaxation of the phase correlation 
(i.e., dephasing) in the basis quantized along an external mag- 
netic field IS [31]. Also, the decoherence of a qubit, which can 
be modeled by a spin- 1/2, is the main obstacle in implement- 
ing scalable quantum computing I4i]. To deal with the spin 
or qubit decoherence, various strategies have been developed, 
including quantum error correction ISllaLZllsl], decoherence- 
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engineering 1251126112711 . In particular, the DD suppresses the 
decoherence by eliminating the qubit-bath coupling through 
stroboscopic rotation of the qubit. An especially interes ting 
DD scheme is the concatenated DD ifisi [igj [20i [III [22i B 
which applies recursively a lower-order periodic pulse se- 
quence as the building block of the next higher order se- 
quence. For an evolution of a short duration T, an Mh or- 
der concatenated DD eliminates the qubit-bath coupling up to 
QrfN+i\ jjjg number of pulses in concatenated DD, how- 
ever, increases exponentially with increasing order A^. Since 
errors are inherently introduced by the controlling pulses, it is 
desirable to have DD sequences with the minimum number of 
controlling pulses. 

An optimal DD scheme was first discovered by Uhrig for a 
pure dephasing spin-boson model 12411 . which uses A^ 7r-pulses 
applied at 



T i = Tsin" 



J^ 



2(N + iy 



for; = l,2,...,N, 



(1) 



to eliminate the dephasing up to Olr'^'^M. Optimal pulse 
sequences for N < 5 have also been noticed by Dhar et 
al. earlier for controlling the Zeno effect \2M- Lee, Witzel 
and Das Sarma conjectured that the Uhrig dynamical decou- 
pling (UDD) may work for a generic pure-dephasing model 
with an analytical verification up to A^ = 9 12^ . Later 
computer-assisted algebra was used to verify the conjecture 
up to A^ = 14 OOll . Aiming at a general proof of the conjec- 
ture, Cardy and Dhar Olll have given a very inspiring though 
unsuccessful attempt by formulating the problem in a time- 



dependent perturbation theory. 

In this Letter, we shall complete the proof of the universal- 
ity of the UDD in suppressing the pure dephasing or the longi- 
tudinal relaxation of a qubit (or spin- 1/2) coupled to a generic 
bath. The proof is based on the observation that to preserve 
the spin coherence up to a given order, one does not have to 
eliminate all terms of the effective qubit-bath coupling to the 
given order as in a generic concatenated DD but just needs to 
eliminate the terms relevant to the decoherence. An extension 
of the proof is that an ideal UDD sequence, for countering the 
longitudinal spin relaxation, can be replaced with a more gen- 
eral UDD sequence consisting of the same ideal Tt-pulses and 
some extra finite-duration pulses as long as certain symmetry 
requirements are fulfilled. 

Ideal UDD for a generic pure dephasing Hamiltonian - Let 
us first consider the ideal UDD pulse sequences for a Hamil- 
tonian of the form 



H^C + &.®Z, 



(2) 



where &. is the qubit Pauli matrix along the z-direction, and 
C and Z are bath operators. This Hamiltonian describes a 
pure dephasing model for it contains no qubit flip processes 
and therefore leads to no longitudinal relaxation but only 
transverse dephasing. Pure dephasing models are of spe- 
cial interest in quantum computing since very often the qubit 
flip terms can be essentially eliminated by applying a strong 
static magnetic field along the z-direction. A specific exam- 
ple is the spin-boson model in which C - YjjUib-bi and 

•Z = Tji(.'^il'^){bi +bi\ with hi being a boson annihilation 
operator It is for this spin-boson model that Uhrig ||24ll has 
discovered the optimal pulse sequences with timing given in 
Eq. dD. 

Now we shall prove that the UDD applies for arbitrary C 
and Z. To overcome the pure dephasing, the ideal UDD se- 
quences consist of 5-pulse 7r-rotations about a transverse axis 
(say, the jc-axis) 12411 . The qubit dephasing is characterized by 
the decay of the expectation value of the raising or lowering 
operator &± = &x ± io'y. 
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where the qubit state-dependent bath propagators under the 



Mh order UDD control are 

j/W ^ ^-i[C±(-\fZ](T-Tft) . . . g-i(C=FZ)(r2-7-,)g-/(C±Z)7-, _ (-4-) 

To show that the Mh order UDD suppresses the pure dephas- 
ing up to O (7"'^+' j for a small T, we just need to prove 



c/i^>'"c/f' = i + o(r'^+'). 



(5) 



By expanding the difference SU = U+ - U- into Taylor series 
and collecting the coefficients term by term, Uhrig has verified 
Eq. ^ for A^ < 14 with computer-assisted algebra JSOll . 

We shall proceed with the formalism of the time-dependent 
perturbation theory due to Cardy and Dhar 13 111 . With the 
standard formulation in the interaction picture, Eq. Q can be 
put in the time-ordered formal expression 



^(N) _ ^-iCT_^^-ij^±FN(t)Z,(t)dt 



(6) 



where £^ is the time-ordering operator, the modulation func- 
tion Ff^it) s (-1)' for t e [Tj,Tj+i] with To = and 
Tfti+i s T, and 
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The difference 5U is given by the Taylor series 

oo 

6U = 2e-'^^ Yu^-if'^' A2,^u (8a) 



i=0 



with 



FN(tn) FN(tn-l)--- ' Ff, (h) 

Jo Jo 

X [Z, (f„ ) Z/ (f„_ 1 ) ■ • ■ Z, (f 1 )] dh dt2--- dt„ . (8b) 

A feature of the expansion of the difference 6U is that it con- 
tains only odd-order terms Aijt+i which are relevant to the de- 
phasing. We just need to show K2k+\ — 0(t^*^\ 
Using the expansion in Eq. (|7]i, we have 

A„ = Yu i^P" ■ ■ ■ Zp.Zp.F,,,,,,...,,T"^'"^P--^''"] , (9) 
[pj] 

where 

is a dimensionless constant independent of T. 
Now the problem is reduced to prove 



F =0 

'^ Pl,P2,-,P„ ^ 



(10) 



for n being odd and « + 'E"i=\ Pj ^ N. For this purpose, we 
make the variable substitution tj - T sin (6*^/2) and define 
the scaled modulation function fN(ff) = FN(Tsirr{6/2)) - 
{-ly for e e [>/(A?-H !),(;■ + l);r/(A^+l)]. With 
sin2''(6i/2) sin = {2i)-^P Yu'JU ^i , sin [(/? - r -n 1) 6*], we can 
write Fp^^p^^...^p^^ as a linear combination of terms of the form 



/. 



,.. ,,„ = de„-\ de. d0, W h [e^ sin (^.-e,) , 

«-' "-^ «-' 'I 



with Ylj\ < Pj + I. Suffices it to show fq^ 



- for odd n 



and 2"=i \qi\ < N. We notice that /n (0) has period of 27t/{N+ 
1) and hence expand it into Fourier series 



//v(0)= Yj ^sm[k(N+l)0]. 



k=l,3,5,- 



(11) 



The key feature of the Fourier expansion to be exploited is 
that it contains only odd harmonics of sin[(A^ + 1)0]. With the 
Fourier expansion, we just need to show that 

d0„ ■ ■ ■ I d02 j d0iY\ COS [rj0j + qj0j) = 0, (12) 

for n being odd, rj being an odd multiple of (A^ -i- 1), and 
2'Li \lj\ ^ ^- With the product-to-sum trigonometric for- 
mula repeatedly used, it can be shown by induction that af- 
ter an even number of variables 0i,02, ■ ■ ■ , 02k have been inte- 
grated over, the resultant integrand as a function of 02k+i can 
be written as the sum of cosine functions of the form 



COS{R2k+\02k+\ + Qlk+lOlk+l), 

with /?2A+i being an odd multiple of (A^ + 1) and \Q2 



(13) 



:k+l 



2k+i\ 



X^2k+ 



qj\. In particular, the last step is 
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cos(R„0„ + QnO„)d0„. 



(14) 



Since R„ is an odd (non-zero, of course) multiple of (N + 1), 
and \Q\ < Y," \qj\ < N, we have R„ + Q„ ^ and the integral 
above must be zero. Thus Eq. ( fTOl ) holds. The proof is done. 

Ideal UDD for suppressing longitudinal spin relaxation - 
Now we consider the most generic qubit-bath Hamiltonian 



H - C + &x ^ X + o-y IS) Y + &r iS> Z, 



(15) 



where o", are the Pauli matrices of the qubit and C, X, Y, and Z 
are bath operators. Without loss of generality, we assume the 
z-axis as the rotation axis for qubit control. We aim to show 
that the spin polarization along the rotation axis Kcf-(r))| is 
preserved up to Ofr'*'^' j under the control of the A^th order 
UDD. The spin polarization is 



\(a,{T)}\ = |<f/W'a-,f/<'^>>| . 



(16) 



where the propagator U can be written as 

jj{N) ^ g-i{C'+(-lfD\(T-TN) . _,^-i{C'-b){T2-Ti)^-i{C'+D)Ti 



(17) 



in which the Hamiltonian has been separated into C' = C + 
a-- (g) Z and b = &,c®X + &y (8> Y. With the definition D/ (f) s 
e''" 'De^''" ', the propagator can be formally expressed as 



U 



(N) 



^-iC'T^^^-if^ FN{l)D,{l)dt 



(18) 



which has the same form as Eq. (|6]l. Following the same 
procedure as for proving Eq. ^, we find that up to O (r'^^' j, 
the expansion of the propagator contains only terms consisting 
of even power of D. Since D contains only the Pauli matrices 
&x and &y and an even power of the two Pauli matrices &"'&/ 
(with itjc + riy being even) is either unity or i&., the propagator 



U 



(N) 



^-iHdTT+OiT"*') 



(19) 



where the effective Hamiltonian i/gft commutes with &-. Thus 
the A^-pulse UDD eliminates the longitudinal qubit relaxation 
uptoO(r^+'). 

UDD with non-ideal pulses: Longitudinal relaxation - 
With the help of Eq. ^^, we realize that Eq. ([TOll holds 
for more general modulation functions F^it) as along as the 
scaled modulation function /Af(0) = F^ \T sin^(0/2)j contains 
only odd harmonics of sin[(A^ + \)ff\ as in Eq. (fTTT l. i.e, 



fN{G) = ^ At sin [{2k + \)i,N + 1)0] . 



(20) 



k=0 



with arbitrary coefficients A^. Motivated by this observation, 
we try to generaUze the UDD to the case of non-ideal pulses. 
Consider the control of the qubit by an arbitrary time- 
dependent magnetic field B{t) applied along the z-direction, 
the general qubit-bath Hamiltonian is 

. 1 
H{t) ^C + &x®X + &y®Y + &;®Z+ -&-B{t). (21) 

In the rotating reference frame following the qubit precession 
under the magnetic field, the Hamiltonian becomes 



HR(t) = C' + cos[(p(t)]D^ + sm[^(t)]D-, 



(22) 



where the precession angle <^(f) - L B{t')dt', C = C+&.®Z, 
D'^ = &x®X+&y®Y, andD" = &x®Y-&y®X. The propagator 
in the rotating reference frame is 



t/ 



-'^'^#exp -/ J Y^^NiODJiDdtl (23) 
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FIG. 1: (color online) An example of (a) the scaled modulation func- 
tions f^(0) for the generalized 3rd order UDD control and (b) the 
coffesponding magnetic field B{t). The dashed lines indicate the cor- 
respondence between the sudden jumps of the modulation function 
f]!/(0) in (a) and the sharp spikes as ideal ;r-pulses in (b). 



the propagator The same way as derive Eq. (|9]l, we find that 
for the «th power of D'^, the expansion in T has coefficients as 



rn+Pi+Pi+--+P„ 



x>-r^rTn^«'wr. 



For n being odd and n + 2'Li Pj ^ N, the multiple integral 
above vanishes [such that the qubit relaxation is suppressed to 
0(t^^^)] as long as the scaled modulation function f^{0) = 

F"^ (t sin"(0/2)j contains only odd harmonics of sin[(A^-i- 1)6] 
as depicted in Eq. ( l20b . This condition is satisfied if and only 
if the scaled modulation functions f^{0) have the following 
symmetries: 

1. periodic with period of2n/{N + l); 

2. anti-symmetric with respect to = J7t/(N + l); 

3. symmetric with respect to = (j + I /2)n/{N + I). 

The anti-symmetry condition requires f^{ff) be either zero or 
discontinuous atO- J7t/{N + 1). But f^{6) and /^(0) cannot 
be simultaneously zero since they have to satisfy the normal- 
ization condition 



[fNi&)f + [fNiO)f 



1, 



(24) 



according to the definition of F^(t). So there must be sud- 
den jumps at least in one of two modulation functions at 
= J7r/(N +1), which means the controlling magnetic field 
with F+(f) = cos[4>{t)], F^{t) = sin[0(f)], and D^(f) = B(0 has to contain a 5-pulse for ;r-rotation at f = Ty. With the 
^iCtj^A^-iC't ^ To consider the qubit relaxation, we just need initial conditions /^(O) - 1 and /^(O) - 0, one can choose 
to examine the odd power of D^ in the Taylor expansion of the field such that ff^{0) is continuous while /^(O) has sudden 



jumps between +1 and -1 at 9 - J7t/(N +1). Thus, a gener- 
alized UDD sequence can be chosen the following way: For 
8 6 [Q,7:/{2N + 2)], f^{6) can be arbitrary but sudden jumps 
from -1 to +1 at = and from +1 to -1 at 7r/{2N + 2), 
and ff^{9) is determined from the normalization condition as 

ff^(0) = + y 1 - [fN{0)f- At other regions, f^(e) are deter- 
mined by the symmetry requirements. The pulse amplitude 
Bit) for the generalized UDD is 

1 d ^ 



Kit) dt 



j=l 



which is a superposition of the ideal UDD pulses and an extra 
component fiextraCO being arbitrary but subject to the symme- 
try requirements. The demand of 5-pulses in the generalized 
UDD is consistent with the previous finding in Ref. 13211 that 
the effect of an ideal 7r-pulse on the evolution of a qubit cou- 
pled to a bath cannot be exactly reproduced by a pulse with a 
finite magnitude. An example of the scaled modulation func- 
tions and the corresponding magnetic field for the generalized 
3rd order UDD control are shown in Fig.[T] Notice that due to 
the variable transformation from 9 to f, the magnetic field Bit) 
does not have the symmetries as the scaled modulation func- 
tions f^i9). For example. Bit) is not periodic and the pulse at 
different time has different width. 

Summary - To summarize, we have proven that with A^ 
ideal 5-pulses for :7r-rotations, the Uhrig dynamical decoupling 
can suppress the pure dephasing or the longitudinal relaxation 
of a qubit (or spin- 1/2) coupled to an arbitrary bath, up to 
QrfN+i\ -pjjg qubit-bath coupling is not eliminated to the 
same order as in generic concatenated dynamical decoupling. 
But the remaining coupling would not induce the qubit deco- 
herence under consideration as it commutes with the relevant 
observable of the qubit. As an extension of the proof, we also 
put forward a design of generalized UDD sequences which 
are the ideal UDD cJ-pulse sequences superimposed with ar- 
bitrary finite-duration pulses satisfying certain symmetry re- 
quirements. It should be pointed out that the present proof of 
the UDD applies either to pure dephasing or longitudinal re- 
laxation and is limited to spin- 1/2. It would be very interesting 
if the UDD can be generalized for simultaneous suppression 
of transverse and longitudinal relaxation and for higher spins. 
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